Naval  Research  Laboratory 

Washington,  DC  20375-5320 


NRL/MR/6390— 15-9665 


Temperature  Histories  of  Structural  Steel 
Laser  and  Hybrid  Laser-GMA  Welds 
Calculated  Using  Multiple  Constraints 


S.G.  Lambrakos 

Center  for  Computational  Materials  Science 
Materials  Science  and  Technology  Division 


December  10,  2015 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway, 
Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 
10-12-2015 

2.  REPORT  TYPE 

NRL  Memorandum  Report 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE  1 

5a.  CONTRACT  NUMBER 

Temperature  Histories  of  Structural  Steel  Laser  and  Hybrid  Laser-GMA  Welds 
Calculated  Using  Multiple  Constraints 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 


5d.  PROJECT  NUMBER 


S.G.  Lambrakos 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory,  Code  6390 
4555  Overlook  Avenue,  SW 
Washington,  DC  20375-5320 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

63-4692-R5521 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


NRL/MR/6390- 15-9665 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR  /  MONITOR’S  ACRONYM(S) 


Office  of  Naval  Research 

One  Liberty  Center 

875  North  Randolph  Street,  Suite  1425 

Arlington,  VA  22203-1995 


ONR 


11 .  SPONSOR  /  MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

Inverse  thermal  analyses  of  structural  steel  deep-penetration  welds  are  presented.  These  analyses  employ  a  methodology  that  is  in  terms  of 
numerical-analytical  basis  functions  and  multiple  constraint  conditions  for  inverse  thermal  analysis  of  steady  state  energy  deposition  in  plate 
structures.  These  analyses  provide  parametric  representations  of  weld  temperature  histories  that  can  be  adopted  as  input  data  to  various  types  of 
computational  procedures,  such  as  those  for  prediction  of  solid-state  phase  transformations  and  mechanical  response.  In  addition,  these  temperature 
histories  can  be  used  to  construct  parametric-function  representations  of  thermal  histories  for  welds  corresponding  to  other  process  parameters 
or  welding  processes  whose  process  conditions  are  within  similar  regimes.  The  present  study  applies  an  inverse  thermal  analysis  procedure 
that  uses  volumetric  constraint  conditions  whose  two-dimensional  projections  are  mappings  onto  transverse  cross  sections  of  experimentally 
measured  solification  boundaries  and  approximate  edges  of  weld  heat-affected  zones. 


15.  SUBJECT  TERMS 

Inverse  Thermal  Analysis  Weld  temperature  histories 

Steel  welds 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Samuel  G.  Lambrakos 

a.  REPORT 

Unclassified 

Unlimited 

b.  ABSTRACT 

Unclassified 

Unlimited 

c.  THIS  PAGE 

Unclassified 

Unlimited 

Unclassified 

Unlimited 

45 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

(202)  767-2601 

1 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.1 8 


Contents 


1.  Introduction . 1 

2.  Inverse  Analysis  Procedure . 2 

3.  Inverse  Thermal  Analysis  of  Structural  Steel  Laser  and  Hybrid  Laser-GMA 

Welds . 3 

4.  Discussion . 4 

5.  Conclusion . 5 

6.  References . 6 


m 


Introduction 


Inverse  thermal  analyses  of  structural  steel  deep-penetration  welds  are  presented.  These 
analyses  provide  a  parameterization  of  temperature  histories  for  prediction  of  properties  within 
the  Heat  Affected  Zone  (HAZ)  of  welds  for  the  regime  considered.  The  present  study  applies  an 
inverse  thermal  analysis  procedure  that  uses  volumetric  constraint  conditions  whose  two- 
dimensional  projections  are  mappings  onto  transverse  cross  sections  of  experimentally  measured 
solidification  and  transformation  boundaries,  or  any  type  of  specified  boundary  within  the 
temperature  field,  including  thermocouple  measurements.  For  the  present  study,  which  considers 
deep-penetration  welds  of  structural  steel,  experimentally  observed  solidification  boundaries  and 
approximate  edges  of  the  heat-affected  zone  (HAZ)  are  adopted  as  temperature-field  constraint 
conditions.  The  inverse  analyses  presented  in  this  study  provide  parametric  representations  of 
weld  temperature  histories  that  can  be  adopted  as  input  data  to  computational  procedures,  such  as 
those  that  predict  solid-state  phase  transformations  and  mechanical  response.  In  addition,  these 
temperature-history  parametric  representations  permit  construction  of  a  model-parameter  space 
that  can  be  adopted  for  inverse  analysis  of  welds  corresponding  to  other  process  parameters  or 
welding  processes  whose  process  conditions  are  within  similar  regimes.  The  construction  of 
temperature  fields  according  to  spatially  and  temporally  distributed  constraint  conditions  using 
linear  combinations  of  basis  functions,  whose  fundamental  trend  features  are  consistent  with 
those  of  the  underlying  process,  represents  a  convenient  approach  to  inverse  analysis  of  energy 
deposition  processes.  Finally,  the  results  of  this  study  contribute  to  a  continually  evolving 
database  of  weld  cross  sections  and  temperature  histories  corresponding  to  specified  weld 
processes,  process  conditions  and  types  of  metals  and  their  alloys. 

The  formal  structure  of  the  inverse  analysis  approach  applied  in  this  study  is  that  of 
parametric  model  representation  of  the  temperature  field  in  terms  of  numerical-analytical  basis 
functions.  The  foundation  of  this  approach  was  introduced  by  [1],  and  provides  a  general 
reduction  of  model  complexity  for  purposes  of  weld  thermal  analysis.  Reduction  of  model 
complexity  is  achieved  by  adopting  numerical-analytical  models  for  either  the  heat  source  or 
temperature  field  (or  both)  in  combination  with  numerical  methods.  Among  notable  studies 
applying  this  approach  are  [2-6].  The  inverse  analysis  approach  applied  here  has  been  used  [7] 
for  construction  of  a  parameter  space  relating  welding  process  parameters,  e.g.,  current  and 
welding  speed,  to  weld  efficiency  for  keyhole  plasma  arc  welding.  Similarly,  the  present  study 
contributes  to  a  parameter  space  relating  shapes  of  weld  solidification  and  transformation 
boundaries,  or  shapes  any  type  of  specified  boundary  within  the  temperature,  to  weld 
temperature  histories,  as  well  as  to  weld  process  parameters.  A  specific  aspect  of  the  inverse 
analysis  methodology  applied  here,  which  permits  convenient  adjustment  of  parameters,  is 
construction  of  discrete  volumetric  source  distributions.  This  aspect  of  the  methodology, 
although  employing  three-dimensional  source  distributions,  is  formally  equivalent  to  the  Myhr- 
Grong  method  [4],  which  typically  employs  discrete  source  distributions  within  two-dimensional 
planes  transverse  to  the  relative  motion  of  heat  source  and  workpiece.  References  [8,9],  and 
reference  therein,  demonstrate  the  use  of  discrete  source  distributions  for  numerical  modeling  of 
transport  phenomena  associated  with  welding  processes. 

The  mathematical  foundation  of  the  inverse  analysis  methodology  for  calculation  of 
temperature  histories  using  linear  combinations  of  numerical-analytical  basis  functions  is  that  of 
least-squares  parameter  optimization  [10-13],  As  discussed  previously  [10-15],  the  inverse 
modeling  approach  compensates  for  lack  of  information  concerning  material  properties  [14-20], 
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This  fundamental  property  of  inverse  models  has  been  discussed  and  demonstrated  previously 
[21-23], 

The  organization  of  the  subject  areas  presented  are  as  follows.  First,  a  brief  description  of 
the  general  procedure  for  inverse  analysis  of  heat  deposition  processes  is  presented.  Second, 
results  of  inverse  thermal  analyses  of  structural  steel  deep  penetration  welds  are  presented.  These 
results  provide  a  quantitative  parametric  representation  of  temperature  histories  for  these  welds 
and  for  any  welds  associated  with  similar  welding  process  conditions.  Third,  a  discussion  is 
presented  concerning  aspects  of  the  inverse  analysis  procedure.  Finally,  a  conclusion  is  given. 

I  nverse  A  nalysis  Procedure 


Following  the  inverse  analysis  approach  [14-20],  a  parametric  model  provides  a  means  for 
the  inclusion  of  information  concerning  the  physical  characteristics  of  a  given  energy  deposition 
process.  Given  the  general  trend  features  of  temperature  fields  associated  with  welding 
processes,  a  consistent  parametric  representation  of  the  temperature  field  for  heat  deposition 
within  structures  characterized  by  a  finite  thickness,  in  terms  of  numerical-analytical  basis 
functions,  is 

Nk  N, 

T(x,t )  =Ta  +^^C(xk)G(x,xk,K,nAt,V)  (Eq  la) 

k=\  n  — 1 


and 


T(xcn,tcn)  =  T: 


where 


G(x,xk,t,K,V )  =  -  exp 


(x-xk-vt)2+(y-yk)2 

4Kt 


(Eq  lb) 
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and  C(x)=^Q(xk)d(x-kk).  (Eq  3) 

k= 1 

where  Q(xk)  is  the  value  of  the  discrete  source  function  at  location  xk.  The  quantities  k,  Vand  l 
are  the  thermal  diffusivity,  welding  speed  and  plate  thickness,  respectively.  The  procedure  for 
inverse  analysis  defined  by  Eqs.(l)-(3)  entails  adjustment  of  the  parameters  C(xk),  xk  and  a? 
defined  over  the  entire  spatial  region  of  the  workpiece.  Formally,  this  procedure  entails 
adjustment  of  the  steady  state  temperature  field  defined  over  the  entire  spatial  region  of  the 
sample  volume.  This  approach  defines  an  optimization  procedure  where  the  temperature  field 
spanning  the  spatial  region  of  the  sample  volume  is  adopted  as  the  quantity  to  be  optimized.  The 
constraint  conditions  are  imposed  on  the  temperature  field  spanning  the  bounded  spatial  domain 
of  the  workpiece  by  minimization  of  the  value  of  the  objective  functions  defined  by 


Zj  = 


(Eq  4) 
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where  t‘  is  the  target  temperature  for  position  xn  =  (. xcn,ycn,zcn ).  The  quantities  wn  ( n=l,...,N )  are 
weight  coefficients  that  specify  relative  levels  of  influence  associated  with  constraint  conditions 

tc  . 

n 

The  inverse  analysis  methodology  defined  by  Eqs.(l)-(4)  is  based  on  a  parametric 
numerical-analytical  model,  which  combines  numerical  integration  with  optimization  of  linear 
combinations  of  numerical-analytical  basis  functions,  which  include  fundamental  solutions  to  the 
heat  conduction  equation  [24]  and  their  Fourier-series  representation.  In  particular,  Eq.(la) 
defines  a  discrete  numerical  integration  over  time,  where  the  time  step  At  is  specified  according 
to  the  average  energy  deposited  during  the  time  At,  for  transition  of  the  temperature  field  to 
steady  state.  It  should  be  noted  that  the  formulation  of  the  inverse  analysis  methodology  defined 
by  Eq.(l)-(4)  is  equipped  with  a  mathematical  structure  that  satisfies  all  boundary  conditions 
associated  with  welding  of  plate  structures  (see  [25]  for  further  discussion). 

In  addition  to  the  parameters  defined  with  respect  to  Eqs.(l)-(4),  the  parameteric  model 
applied  for  inverse  analysis  includes  a  length  scale  parameter  ls,  where  in  general  ls  <  l  defined 
by  Eq.(2),  for  specification  of  the  spatial  scale  of  the  calculated  temperature  field  with  respect  to 
which  parameters  are  adjusted.  This  length  scale  parameter  provides  for  inclusion  of  more  details 
of  shape  features  of  measured  solidification  boundaries  to  be  adopted  as  constraint  conditions. 
However,  for  the  analyses  that  follow,  the  two  length  parameters,  i.e.,  the  plate  thickness  l  and 
the  depth  ls  of  the  speciifed  region  of  the  temperature  field  to  be  calculated,  are  taken  to  be  the 
same. 

Inverse  Thermal  Analysis  of  Structural  Steel  Laser  and  Hybrid  Laser-GMAWelds 

In  this  section  results  of  inverse  thermal  analyses  of  structural  steel  deep-penetration  welds 
are  described,  which  correspond  to  different  weld  process  conditions  and  associated  process- 
control  parameters.  The  significance  of  the  inverse-problem  approach  for  these  analyses  is  that 
the  nature  of  the  coupling  of  the  energy  source  to  the  workpiece,  which  is  a  function  of  beam 
power  and  process  control  parameters,  is  in  principle  difficult  to  specify  relative  to  analysis 
based  on  the  direct-problem  approach. 

Previous  studies  considered  inverse  thermal  analysis  of  steel  welds  using  single  cross  section 
[25]  and  volumetric  [26]  constraint  conditions.  A  series  of  studies  [27-29]  considered  inverse 
thermal  analysis  of  welds  using  multiple  constraint  conditions,  where  both  solidification  and 
phase-transformation  boundaries  were  adopted  as  constraint  conditions.  In  particular,  welds  of 
CP-Ti  and  a  Ti-6A I-4V  were  considered  [30-32],  where  weld-cross-section  measurements  of 
both  solidifcation  and  estimated  a-/3  phase-transformation  boundaries  were  adopted  as 
constraint  conditions  on  the  calculated  temperature  field.  In  these  studies,  motivation  for 
adopting  a-/3 phase  transformation  boundaries  as  constraint  conditions  was  that  in  practice,  for 
welds  of  Ti  and  its  alloys,  one  can  associate  (approximately)  this  boundary  with  the  observed 
edge  of  the  HAZ,  and  accordingly,  specify  an  isothermal  boundary  of  known  temperature.  The 
present  study  is  based  on  a  similar  motivation.  In  particular,  for  steel  welds,  especially  those  of 
structural  steel,  the  HAZ  consists  by  various  regions  whose  microstructures  have  been  well 
characterized  (see  reference  [33]  for  further  discussion).  These  regions  are  described 
schematically  in  Fig.  1.  Referring  to  this  figure,  it  should  be  noted  that,  for  purposes  of  inverse 
thermal  analysis,  one  should  be  able  to  associated  the  observed  edge  of  the  HAZ  (approximately) 
with  a  maximum  temperature  700  °C.  As  discussed  in  reference  [27],  for  regions  within  the 
workpiece  not  close  to  the  energy  source,  reasonable  estimates  of  temprature-field  values,  at 
specified  locations,  are  sufficient  for  imposing  constraint  conditions. 
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The  structural  steel  deep-penetration  welds,  whose  inverse  analysis  is  presented  here, 
consist  of  laser  and  laser-GMA  hybrid  welds  [34],  The  analyses  presented  here  entail  calculation 
of  the  steady  state  temperature  field  for  a  specified  range  of  sizes  and  shapes  of  inner  surface 
boundaries  s,  defined  by  the  solidification  boundary  and  experimentally  observed  estimate  of  the 
HAZ  edge.  The  shapes  of  these  boundaries  are  determined  experimentally  by  analysis  of 
transverse  weld  cross  sections  showing  microstructure  revealing  solidification  and  estimated 
HAZ-edge  boundaries.  For  calculations  adopting  the  solidification  boundaries  as  constraints,  the 
parameter  values  assumed  are  k=  5.88  x  1CT6  m2s'\  TM  =  1503.0  °C.  For  calculations  adopting 
the  HAZ-edge  boundaries  as  constraints,  the  parameter  values  assumed  are  k  =  5.88  x  10'6  mV1, 
The  =  700  °C.  As  discussed  previously  [25,  27],  reasonable  estimates  of  k  and  Tu  are  sufficient 
for  inverse  analysis.  For  this  study,  the  same  thermal  diffusivity  is  assumed  for  calculation  of 
temperature  fields  using  solidification  and  estimated  HAZ-edge  constraints.  This  assumption  is 
sufficient,  within  reasonable  estimates,  in  that  the  set  of  parameters  C(xk),  k=l,...,Nk,  and k are 
not  uniquely  determined  by  inverse  analysis.  Thus,  changing  estimated  values  of  k  would 
require  different  values  of  c(xk)  in  order  to  satisfy  specified  constraint  conditions  associated 
with  tm  and  THEl.  With  respect  to  inverse  analysis,  the  interpretation  of  k  as  both  an  estimated 
material  property  and  adjustable  parameter  is  emphasized  within  the  following. 

The  goal  of  the  present  analysis  is  determination  of  a  set  of  parameters  that  can  serve  as 
initial  estimates  for  parameter  adjustment  with  respect  to  deep  penetration  welds  of  steels,  whose 
process  parameters  are  within  similar  regimes.  Parameter  adjustment  with  respect  to  other  welds, 
which  assume  the  results  of  this  study  as  initial  estimates,  would  adopt  k  and  tm  as  adjustable 
parameters,  as  well  as  the  discrete  source  function  C{xk).  Values  of  the  workpiece  thickness  l 
and  welding  speed  V  defined  in  Eq.  (2)  are  given  in  the  figures  below.  The  upstream  boundary 
constraints  on  the  temperature  field,  Tc=  Tm  and  Tc=  THe  for  (yc,zc)  defined  in  Eq.  (lb),  are 
given  in  Table  1  for  the  solidification  and  HAZ-edge  boundaries,  respectively.  Given  in  Tables  2 
through  17  are  values  of  the  discrete  source  function  that  have  been  calculated  according  to  the 
constraint  conditions  and  weld  process  specifications  given  in  Table  1.  The  relative  location  of 
each  discrete  source  is  specified  according  to  Fig.  2.  Shown  in  Figs.  3  through  42  are 
experimentally  measured  transverse  weld  cross  sections  of  solidification  and  HAZ-edge 
boundaries  [34],  and  different  planer  slices  of  the  steady  state  temperature  field  that  have  been 
calculated  according  to  the  constraint  conditions  given  in  Table  1.  Referring  to  the  planer  slices 
of  the  calculated  temperature  fields  shown  in  these  figures,  it  can  be  seen  that  all  boundary 
conditions  are  satisfied,  namely  the  condition  T(x,t)=TMor  T(x,t)  =  THE  at  the  solidification  or 
HAZ-edge  boundary,  respectively,  and  VT-n  =  o  at  surface  boundaries,  where  n  is  normal  to  the 
surface. 

Discussion 

The  inverse  analysis  procedure  entails  calculating  either  a  three-dimensional  solidification  or 
HAZ-edge  boundary  using  experimentally  measured  constraint  conditions,  and  the  temperature 
field  consistent  with  the  isothermal  surface  associated  with  this  boundary. 

Shown  in  Figs.  4,  9,  14,  19,  24,  29,  34  and  39  are  two-dimensional  slices  of  the  calculated 
three-dimensional  temperature  field  obtained  using  the  constraint  conditions  given  in  T able  1  for 
the  measured  solidification  boundary,  which  are  parallel  to  the  relative  motion  of  laser  or  laser- 
GMA  source  and  workpiece.  Shown  in  Figs.  5, 10, 15,  20,  25,  30,  35  and  40  are  two-dimensional 
slices  of  this  three-dimensional  temperature  field  that  are  perpendicular  to  the  relative  motion  of 
laser  or  laser-GMA  source  and  workpiece.  Referring  to  these  figures,  it  should  be  noted  that  ^  = 
0  has  been  assigned  arbitrarily  to  a  two-dimensional  slice  at  the  leading  edge  of  the  solidification 
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boundary.  Accordingly,  shown  in  Figs.  5,  10,  15,  20,  25,  30,  35  and  40  is  passage  with  time  of 
the  calculated  three-dimensional  solidification  boundaries  through  experimentally  measured 
transverse  cross  sections  of  these  boundaries. 

Shown  in  Figs.  6, 11, 16,  21,  26,  31,  36  and  41  are  two-dimensional  slices  of  the  calculated 
three-dimensional  temperature  field  obtained  using  the  constraint  conditions  given  in  Tables  1 
through  4  for  the  estimated  HAZ-edge  boundary,  which  are  parallel  to  the  relative  motion  of 
laser  or  laser-GM  A  source  and  workpiece.  Shown  in  Figs.  7,  12,  17,  22,  27,  32,  37  and  42  are 
two-dimensional  slices  of  this  three-dimensional  temperature  field  that  are  perpendicular  to  the 
relative  motion  of  laser  or  laser-GM  A  source  and  workpiece.  Again,  referring  to  these  figures,  it 
should  be  noted  that  t  =  0  has  been  assigned  arbitrarily  to  a  two-dimensional  slice  at  the  leading 
edge  of  the  solidification  boundary.  Again,  accordingly,  shown  in  Figs.  7,  12, 17,  22,  27,  32,  37 
and  42  is  passage  with  time  of  the  calculated  three-dimensional  solidification  boundaries  through 
experimentally  measured  transverse  cross  sections  of  these  boundaries. 

Referring  to  the  calculated  temperature  fields  shown  above,  it  is  important  to  note  that  the 
constraint  conditions  on  the  calculated  three-dimensional  solidification  and  HAZ-edge 
boundaries  are  that  the  projections  of  all  their  two-dimenional  transverse  slices,  as  a  function  of 
time,  are  consistent  with  the  experimentally  measured  transverse  cross  sections  of  the 
solidification  and  transformation  boundaries,  respectively. 

The  final  stage  of  the  analysis  entails  calculation  of  temperature  histories  as  a  function  of 
transverse  position  within  the  cross  section  of  the  weld.  This  is  accomplished  by  constructing  a 
steady  state  three-dimensional  temperature  field  that  is  consistent  with  experimentally  measured 
constraint  conditions  for  the  solidification  boundary  or  for  both  solidification  and  HAZ-edge 
boundaries.  In  cases  where  both  the  solidification  and  HAZ-edge  boundaries  are  adopted  as 
constraints,  this  construction  combines  two  regions  of  the  temperature  field.  One  region  consists 
of  the  calculated  temperature  field  obtained  using  HAZ-edge  boundary  constraints  for  all 
temperatures  less  than  and  equal  to  the  temperature  tHe ■  The  other  region  consists  of  the 
temperature  field  that  is  obtained  by  three-dimensional  interpolation  between  isothermal  surfaces 
tm  and  tHe,  which  are  associated  with  the  solidification  and  HAZ-edge  boundaries,  respectively. 
In  general,  the  numerical  procedure  for  interpolation  between  constrained  isotherms  is  a  separate 
issue  for  consideration. 

Conclusion 

A  specific  objective  of  this  report  is  to  further  examine,  for  the  case  of  structural  steel  welds, 
the  concept  of  using  HAZ-edge  boundary  constraints  for  inverse  thermal  analysis.  This  concept 
was  introduced  previously  for  inverse  thermal  analysis  of  CP-Ti  and  a  Ti-6A  I-4V  welds  [30-32]. 
A  general  objective  of  this  report  is  to  describe  a  quantitative  inverse  thermal  analysis  of 
structural  steel  deep-penetration  welds  corresponding  to  various  weld  process  parameters  and  to 
construct  numerical-analytical  basis  functions  that  can  be  used  by  weld  analyst  to  calculated 
weld  temperature  histories,  which  are  for  welding  processes  associated  with  similar  process 
conditions.  This  report  contributes  to  the  continuing  evolution  of  a  parametric  representation  of 
the  temperature  field  T(x,t,K,v,  Si,  l)  for  inverse  thermal  analysis  of  welds  associated  with 
different  types  of  metals,  their  alloys  and  weld  process  conditions.  The  weld  temperature 
histories  obtained  by  inverse  analysis  could  in  practice  be  used  to  predict  solid-  state  phase 
transitions  and  mechanical  response  characteristics  associated  with  temperature  histories  within 
localized  spatial  regions,  and  thus  with  evolution  of  plastic  and  elastic  strains,  resulting  in 
distortion  and  residual  stresses. 
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Table  1  Temperature  field  constraint  conditions  at  positions  (y«zc)  at  solidification  and  HAZ- 
edge  boundaries  on  transverse  cross  sections  of  structural  steel  laser  and  hybrid  laser-GTA  welds 

Solidification  Boundaries 


Solidification  Boundaries 


HAZ  Edge  Boundaries 


HAZ  Edge  Boundaries 


WELD  5 

WELD  6 

WELD  7 

WELD  8 

(10.18,  0.2679) 

(10.47,  0.2830) 

(12.32,  0.3159) 

(12.38,  0.3175) 

(9.912,  1.340) 

(10.47,  1.415) 

(11.69,  1.580) 

(12.38,3.175) 

(9.644,  2.679) 

(9.905,  2.830) 

(10.11, 3.159) 

(11.75,4.763) 

(8.573,4.019) 

(9.056,  4.019) 

(8.845,  4.739) 

(10.80,  6.35) 

(7.501, 5.358) 

(7.924,  5.66) 

(6.95,  6.318) 

(9.842,  7.938) 

(5.894,  6.698) 

(6.509,  7.075) 

(5.686,  7.898) 

(8.255,  9.525) 

(4.554,  8.037) 

(5.094,  8.037) 

(4.423,  9.477) 

(6.35,  12.7) 

(4.019,  9.377) 

(3.962,  1 1 .32) 

(4.107,  11.06) 

(5.398,  15.89) 

(3.751,  10.72) 

(3.962,  12.74) 

(3.475,  12.64) 

(5.715,  19.05) 

(3.483,  13.40) 

(3.962,  14.15) 

(3.159,  14.22) 

(6.033,  20.0) 

(3.483,  14.73) 

(3.962,  15.0) 

(3.475,  15.80) 

(3.791,  17.37) 

(3.791,  18.0) 

Table  2  Volumetric  source  function  c(xk)  calculated  according  to  solidification-boundary 
constraint  conditions  given  in  Table  1,  where  M  =  (12/60)  mm,  xk  =  yk  =  0.0  for  k  =  1  to  13 
(WELD  1). 
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8 

0.07 
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9 
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CD 

i — 1 
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45 

11 

0.06 

50 

12 

0.05 

55 

13 

0.05 

60 

Table  3  Volumetric  source  function  c(xk)  calculated  according  to  HAZ-edge  constraint 
conditions  given  in  T able  1,  where  m  =(12/60)  mm,x^  =  ^  =  0.0for^  =  lto  13  (WELD  1). 
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0.04 

55 

13 
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60 

Table  4  Volumetric  source  function  c(xk)  calculated  according  to  solidification-boundary 
constraint  conditions  given  in  Table  1,  where  M  =  (12/60)  mm,  xk  =  yk  =  0.0  for  k  =  1  to  13  and 
z*  =  1  for  k  =  14  to  17  (W  E  L  D  2). 
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0.09 
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xk 


M  ) 


y±  (  A/) 


14 
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-10.0 


0.0 
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2.4 


10.0 


0.0 
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2.4 


0.0 
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17 


2.4 


0.0 
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Table  5  Volumetric  source  function  c(xk)  calculated  according  to  HAZ-edge  constraint 
conditions  given  in  Table  1,  where  M  =  (12/60)  mm,  xk  =  yk  =  0.0  for  k  =  1  to  13  and  zk  =  1  for  k 
=  14  to  17  (WELD  2). 
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%k  (  ) 

W  (  A/) 
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-10.0 
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10.0 

0.0 
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1.5 

0.0 
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1.5 

0.0 

10.0 

Table  6  Volumetric  source  function  c(xk)  calculated  according  to  solidification-boundary 
constraint  conditions  given  in  Table  1,  where  m  =  (15/60)  mm,  xk  =  yk  =  0.0  for  k  =  1  to  13 
(WELD  3). 
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60 

Table  7  Volumetric  source  function  c(xk)  calculated  according  to  HAZ-edge  constraint 
conditions  given  in  T able  1,  where  a/  =(15/60)  mm,;c*  =  ;y*  =  0.0for&  =  lto  13  (WELD  3). 
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Table  8  Volumetric  source  function  C(xk)  calculated  according  to  solidification-boundary 
constraint  conditions  given  in  Table  1,  where  M  =  (15/60)  mm,  xk  =  yk  =  0.0  for  k  =  1  to  13 
(WELD  4). 
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Table  9  Volumetric  source  function  c(xk)  calculated  according  to  HAZ-edge  constraint 
conditions  given  in  T able  1,  where  m  =(15/60)  mm,x^  =  ^  =  0.0for^  =  lto  13  (WELD  4). 
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Table  10  Volumetric  source  function  c(xk)  calculated  according  to  solidification-boundary 
constraint  conditions  given  in  Table  1,  where  M  =  (15/60)  mm,  xk  =  yk  =  0.0  for  k  =  1  to  13 
(WELD  5). 
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60 

Table  11  Volumetric  source  function  c(xk)  calculated  according  to  HAZ-edge  constraint 
conditions  given  in  T able  1,  where  m  =(15/60)  mm,x^  =  ^  =  0.0for^  =  lto  13  (WELD  5). 
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60 

Table  12  Volumetric  source  function  c(xk)  calculated  according  to  solidification-boundary 
constraint  conditions  given  in  Table  1,  where  M  =  (3.0/60)  mm,  xk  =  yk  =  0.0  for  k  =  1  to  13 
(WELD  6). 
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Table  13  Volumetric  source  function  c(xk)  calculated  according  to  HAZ-edge  constraint 
conditions  given  in  T able  1,  where  m  =(15/60)  mm,x^  =  ^  =  0.0for^  =  lto  13  (WELD  6). 
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Table  14  Volumetric  source  function  c(xk)  calculated  according  to  solidification-boundary 
constraint  conditions  given  in  Table  1,  where  m  =  (18/60)  mm,  xk  =  yk  =  0.0  for  k  =  1  to  13 
(WELD  7). 
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Table  15  Volumetric  source  function  c(xk )  calculated  according  to  HAZ-edge  constraint 
conditions  given  in  T able  1,  where  M  =(18/60)  mm,x^  =  ^  =  0.0for^  =  lto  13  (WELD  7). 
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Table  16  Volumetric  source  function  c(xk)  calculated  according  to  solidification-boundary 
constraint  conditions  given  in  Table  1,  where  M  =  (20/60)  mm,  xk  =  yk  =  0.0  for  k  =  1  to  13 
(WELD  8). 


k 

cy  2  5.0 

z*  (A/) 

i 

0.40 

1 

2 

0.38 

5 

3 

0.35 

10 

4 

0.15 

15 

5 

0.15 

20 

6 

0.15 

25 

7 

0.15 

30 

8 

0.14 

35 

9 

0.12 

40 

CD 

i — 1 

0.12 

45 

11 

0.12 

50 

12 

0.12 

55 

13 

0.09 

60 

Table  17  Volumetric  source  function  c(xk)  calculated  according  to  HAZ-edge  constraint 
conditions  given  in  T able  1,  where  m  =(20/60)  mm,x^  =  ^  =  0.0for^  =  lto  13  (WELD  8). 
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Figure  1.  Schematic  representation  of  different  regions  within  the  HAZ  for  steel  welds  and 
approximate  location  of  experimentally  observable  HAZ  edge. 


Figure  2.  Indexing  scheme  for  relative  locations  of  discrete  sources  C(xk) ,  k=l,...,Nk. 
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Figure  3.  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and  HAZ- 
edge  boundaries  for  steel  laser  weld,  where  V  =1.333  cm/s  (W  eld  1). 
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Figure  4.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  solidification  boundary  (Weld  1). 
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Figure  5.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  solidification  boundary,  where  a?=az/V,  m  =  (12/60) 
mm  and  V  =  1.333  cm/s  (W eld  1). 
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Figure  6.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  HAZ-edge  boundary  (Weld  1). 
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Figure  7.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  HAZ-edge  boundary,  where  a?=a//V,  m  =  (12/60) 
mm  and  V  =  1.333  cm/s  (W eld  1). 
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Figure  8.  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and 
transformation  (HAZ-edge)  boundaries  for  steel  laser-GMA  hybrid  weld,  where  V=2. 333  cm/s 
(Weld  2). 
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Figure  9.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section  at 
symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  solidification  boundary  (Weld  2). 
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Figure  10.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  solidification  boundary,  where  a?=az/V,  m  =  (12/60) 
mm  and  V  =  2.333  cm/s  (W eld  2). 
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Figure  11.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  HAZ-edge  boundary  (Weld  2). 


23 


i  §  i  i  §  !  2 

Figure  12.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  HAZ-edge  boundary,  where  At  =a//V,  M  =  (12/60) 
mm  and  V  =  2.333  cm/s  (W eld  2). 
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Figure  13.  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and 
transformation  (HAZ-edge)  boundaries  for  steel  laser-GMA  hybrid  weld,  where  V  =1.8333  cm/s 
(Weld  3). 


Figure  14.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  solidification  boundary  (Weld  3). 


25 


t  =  20A  t 


t  =  25A  t 


T{°C) 


Figure  15.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  solidification  boundary,  where  m  =mN ,  m  =  (15/60) 
mm  and  V  =  1.8333  cm/s  (W eld  3). 


Figure  16.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  HAZ-edge  boundary  (Weld  3). 
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Figure  17.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  HAZ-edge  boundary,  where  a?=az/V,  m  =  (15/60) 
mm  and  V  =  1.8333  cm/s  (W eld  3). 
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Figure  18.  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and 
transformation  (HAZ-edge)  boundaries  for  steel  laser-GMA  hybrid  weld,  where  V  =2.1666  cm/s 
(Weld  4). 
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Figure  19.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  solidification  boundary  (Weld  4). 
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Figure  20.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  solidification  boundary,  where  m  =a//V ,  m  =  (15/60) 
mm  and  V  =  2.1666  cm/s  (W eld  4). 
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Figure  21.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  HAZ-edge  boundary  (Weld  4). 
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Figure  22.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  HAZ-edge  boundary,  where  m  =a//V ,  M  =  (15/60) 
mm  and  V  =  2.1666  cm/s  (W eld  4). 
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Figure  23.  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and 
transformation  (HAZ-edge)  boundaries  for  steel  laser-GMA  hybrid  weld,  where  V=1.5  cm/s 
(Weld  5). 
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Figure  24.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  solidification  boundary  (Weld  5). 
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Figure  25.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  solidification  boundary,  where  m  =a//V ,  M  =  (15/60) 
mm  and  V  =  1.5  cm/s  (W eld  5). 


Figure  26.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  HAZ-edge  boundary  (Weld  5). 
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Figure  27.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  HAZ-edge  boundary,  where  a?=a//V,  a/  =  (15/60) 
mm  and  V  =  1.5  cm/s  (W eld  5). 
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Figure  28.  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and 
transformation  (HAZ-edge)  boundaries  for  steel  laser-GMA  hybrid  weld,  where  V=1.5  cm/s 
(Weld  6). 
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Figure  29.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  solidification  boundary  (Weld  6). 
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Figure  30.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  solidification  boundary,  where  m  =ai/V ,  m  =  (15/60) 
mm  and  V  =  1.5  cm/s  (W eld  6). 


Figure  31.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  HAZ-edge  boundary  (Weld  6). 
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Figure  32.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  HAZ-edge  boundary,  where  a?=a//V,  m  =  (15/60) 
mm  and  V  =  1.5  cm/s  (W eld  6). 
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Figure  33.  Experimentally  measured  transverse  weld  cross  sections  of  solidification  and 
transformation  (HAZ-edge)  boundaries  for  steel  laser-GMA  hybrid  weld,  where  V=1.5  cm/s 
(Weld  7). 
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Figure  34.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  solidification  boundary  (Weld  7). 


t  =  0  t  =  lOAt 


t  =  20  At  t  =  40  At 


8  8  8  5  |  g  | 


Figure  35.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  solidification  boundary,  where  a?=a//V,  m  =  (18/60) 
mm  and  V  =  1.5  cm/s  (W eld  7). 
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Figure  36.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  HAZ-edge  boundary  (Weld  7). 
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Figure  37.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  HAZ-edge  boundary,  where  a?=a//V,  m  =  (20/60) 
mm  and  V  =  1.5  cm/s  (W eld  7). 
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(Weld  8). 


Figure  39.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  solidification  boundary  (Weld  8). 
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Figure  40.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  solidification  boundary,  where  a?=az/V,  m  =  (20/60) 
mm  andV  =8.333  mm/s  (Weld  8). 
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Figure  41.  Two-dimensional  slices,  at  half  workpiece  top  surface  and  longitudinal  cross  section 
at  symmetry  plane,  of  three-dimensional  temperature  field  (°C)  calculated  using  cross  section 
information  given  in  Table  1  for  HAZ-edge  boundary  (Weld  8). 
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Figure  42.  Temperature  history  (°C)  of  transverse  cross  section  of  weld  calculated  using  cross 
section  information  given  in  Table  1  for  HAZ-edge  boundary,  where  a?=a//V,  m  =  (20/60) 
mm  andV  =8.333  mm/s  (Weld  8). 
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